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Center,  Purdue  University,  West  Lafayette,  Indiana,  under  USAF  Contract 
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to  30  June  1970.  The  contract  was  administered  by  the  Materials  Physics  Division, 
Air  Force  Materials  Laboratory  Wright  Patterson  Air  Force  Base,  Ohio,  with 
Mr.  Hyman  Marcus  as  project  engineer.  The  report  was  submitted  by  the  author 
for  publication  as  an  AFML  technical  report  on  20  July  1970. 
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ABSTRACT 


The  development  of  the  multiproperty  apparatus  to  include  high  temperature 
determinations  of  specific  heat,  enthalpy  and  thermal  diffusivity  in  addition  to  the 
properties  already  measured  is  described.  The  latter  properties  include  thermal 
conductivity,  electrical  resistivity,  Lorenz  ratio,  total  hemispherical  emittance, 
spectral  emittance  at  0.65  microns  and  Thomson  coefficient.  The  specific  heat  de¬ 
termination  involves  small  transients  from  an  established  steady-state  temperature. 
The  accuracy  and  reproducibility  of  the  specific  heat  determination  is  controlled  by 
the  accuracy  obtained  in  measuring  the  rate  of  change  of  temperature.  When  photo¬ 
graphing  an  oscilloscopic  trace  of  an  automatic  pyrometer  output,  the  inaccuracy 
was  about  8  percent.  By  using  a  digital  data  acquisition  system,  it  is  believed  that 
data  accurate  within  1  percent  can  be  achieved  at  high  temperatures. 
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NOMENCLATURE 


Definition 

Cross-sectional  area  perpendicular  to  Z-axis 
Brightness  temperature 
Specific  heat 
Electric  current 

Electric  current  through  infinitely  long  sample 
Joules 

Circumference  of  sample 

Temperature  of  the  sample  at  any  location  Z 

Rate  of  change  of  temperature 

True  temperature 

Ambient  temperature 

Temperature  at  center  of  infinitely  long  sample 

Voltage  drop 

Watts 

Length  coordinate  along  longitudinal  axis  of  sample 

Total  hemispherical  absorptance 
Total  hemispherical  emittance 
Spectral  emittance  at  0.  65  microns 
Density 

Thermal  conductivity 
Thomson  coefficient 
Electrical  resistivity 
Stefan-Boltzmann  constant 


I.  INTRODUCTION 


There  are  two  common  methods  of  measuring  the  thermal  conductivity  of 
solids  at  high  temperatures.  These  are  the  radial  heat  flow  method  in  which  a 
calibrated  heater  or  heat  sink  is  inserted  in  the  center  of  a  tubular  sample  and  the 
radial  temperature  gradient  determined,  and  the  direct  heating  method  in  which  a 
long  thin  sample  is  clamped  between  two  electrodes  and  heated  by  passing  'egulated 
current  through  it  while  determining  the  longitudinal  temperature  gradient.  The  di¬ 
rect  heating  method  is  simpler  experimentally  but  more  complex  mathematically. 


For  example  the  differential  equation  describing  energy  transport  in  a  long 
thin  rod  at  steady  state  heated  electrically  in  vacuo  and  free  to  radiate  to  its 
enclosure  is 


jL,' 

dZ  \ 


+  lie  _  (T<  _  T 

A2  A  U  o 


I  dT  _  . 
”  4  A  dZ  0 


Eq.  (1) 


where  I  is  the  DC  current,  A  is  the  cross-sectional  area,  P  is  the  circumference, 
CT  is  the  Stefan-Britzmann  constant,  T  is  the  temperature  at  any  longitudinal  posi¬ 
tion  Z,  Tq  is  the  temperature  of  the  enclosure,  A  is  the  thermal  conductivity,  p  is 
the  electrical  resistivity,  is  the  total  hemispherical  emittance  and  p  is  the 
Thomson  coefficient.  There  are  several  important  things  to  note  in  Eq.  (1).  First, 
X,  p,  €pj,  and  y  are  functions  of  temperature  (and  therefore  of  position).  There  is 
no  known  closed  solution  for  X  even  for  the  relatively  simple  case  in  which  u  =  0 
and  the  temperature  dependencies  of  X,p  and  6^  are  negligible.  The  experimentalist 
measures  T  as  a  function  of  Z,  whereas  Eq.  (1)  contains  dT/dZ  and  d2T/dZ2  in 
addition  to  powers  of  T.  Therefore,  it  is  necessary  to  employ  very  restrictive 
mathematical  approximations  in  order  to  linearize  Eq.  (1),  use  numerical  integra¬ 
tion  techniques  or  differentiate  an  appropriate  expression  for  T(Z)  to  obtain  dT/dZ 
and  d2T/dZ2  in  order  to  obtain  X  from  the  experimental  data.  The  various  variants 
have  been  described  in  several  papers  Taylor  and  coworkers  1  *  ^ 

have  developed  several  analytical  schemes  which  permit  the  calculation  of  X  from 
the  experimental  data  without  the  use  of  restrictive  mathematical  approximations. 

Using  these  techniques  they  have  obtained  very  accurate  data  for  X,  o  and  £„  to 

(2  5  8  91  " 

high  temperatures  ’  ’  ’ '  \  For  example,  portions  of  the  p  and  £„  curves  for 

ir))  H 

arc-cast  tungsten  are  shown  in  Figures  1  and  2  while  the  X  results  '  are 

shown  in  Figure  .3. 

I, 


Taylor  et  al.  ^  showed  the  necessity  of  measuring  p,  and  X  on  the 

same  sample  in  the  same  apparatus.  Basically  the  reason  for  this  is  that  small 

uncertainties  in  p  or  €„  cause  a  large  uncertainty  in  X  because  the  second  and 

n 

third  terms  of  Eq.  (1)  are  of  nearly  equal  magnitude  and  the  fourth  term  is  rela¬ 
tively  small.  Thus  the  magnitude  of  the  term  containing  X  is  essentially  the  dif¬ 
ference  between  the  magnitudes  of  nearly  equal  terms.  This  is  illustrated  in 
Table  I  in  which  the  magnitudes  of  the  various  terms  at  selected  positions  along  a 
typical  profile  are  computed.  At  the  position  Z  =  22. 642  cm,  it  is  seen  that  the 
Joulean  heating  term  I2p/A2  is  677. 165  watts/cm8,  the  radiation  term  is  640.662 
watts/cm8  and  the  Thomson  term  is  1.021  watts/cm8.  The  term  containing  X  is 
only  35.481  watts/cm8  and  thus  a  small  error  in  p  or  has  a  large  effect  on  this 
term.  Because  of  this  situation,  the  successful  employment  of  the  direct  heating 
method  involves  multiproperty  determinations. 

While  the  muuiematical  techniques  have  not  been  optimized  as  yet  for 
determining  p,  it  appears  feasible  to  do  so  and  Taylor  and  coworkers  are  pursuing 
this  objective  Also  the  spectral  emittance  at  0.65  microns  is  readily  deter¬ 
mined  from  the  experimental  data  generated  during  a  direct  heating  experiment. 
Since  X  and  p  are  determined  on  the  same  specimen,  the  Lorenz  function  can  also 
be  calculated  from  the  data.  Thus  the  direct  heating  method  has  been  used  to 
measure  six  properties  simultaneously  on  the  same  sample. 

Because  of  the  ideal  sample  geometry,  it  appears  feasible  to  measure  at 
least  seven  additional  thermophysical  properties  on  the  same  sample  with  accuracies 
comparable  to  those  obtainable  under  optimum  conditions  for  each  property.  These 
additional  properties  are  thermal  expansion,  specific  heat,  enthalpy,  thermal  dif- 
fusivity,  Seebeck  coefficient,  Peltier  coefficient,  and  Richardson  coefficient.  The 
present  project  has  as  its  goal  the  extension  of  the  direct  heating  technique  to  in¬ 
clude  specific  heat  determinations.  Since  the  enthalpy  can  be  obtained  by  merely 
integrating  the  specific  heat  results,  it  would  be  calculable  fron  the  specific  heat 
data.  Similarly,  because  the  thermal  conductivity  would  be  determined  on  the 
same  specimen,  measuring  the  specific  heat  would  permit  the  calculation  of  thermal 
diffusivity,  provided  the  density  is  known. 
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TEMPERATURE  (K) 

FIGURE  I.  ELECTRICAL  RESISTIVITY  OF  TUNGSTEN 
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TEMPERATURE,  K 

FIGURE  2.  TOTAL  HEMISPHERICAL  EMMITTANCE  OF  TUNGSTEN 
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FIGURE  3.  THERMAL  CONDUCTIVITY  OF  TUNGSTEN 


MAGNITUDES  OF  CONDUCTION,  HEATING,  RADIATION,  AND  THOMSON 
TERMS  AT  SELECTED  POSITIONS  ALONG  A  TYPICAL  TEMPERATURE  PROFILE 
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II.  BACKGROUND  ON  SPECIFIC  HEAT  DETERMINATIONS 


The  direct  determination  of  specific  heat  (Cp)  at  high  temperatures  has  been 
made  using  pulse  heating  or  modulating  current  techniques  14^.  The 

specific  heat  has  also  been  determined  at  high  temperatures  by  measuring  the 
enthalpy  as  a  function  of  temperature  by  drop -calorimetry  and  differentiating.  The 
results  of  these  various  approaches  are  often  in  good  agreement  from  room  tem¬ 
perature  to  temperatures  about  2/3  of  the  melting  point  but  often  diverge  markedly 
as  the  temperature  approaches  the  melting  point  (10*  For  example,  compare 
the  data  of  Hoch  and  Johnston  (drop  calorimetry)  and  Taylor  and  Finch 
(pulse  heating)  for  tantalum  given  in  Table  II.  The  results  obtained  by  the  modula- 
tion  method  by  Lowenthal '  '  agree  with  the  Taylor  and  Finch  results,  whereas 

the  modulation  results  of  Kraftmahker  agree  with  Hoch  and  Johnston. 


TABLE  II 

SPECIFIC  HEAT  OF  TANTALUM  (J/gm  K) 


T  (K) 

Hoch  and 
Johnstont16) 

Taylor  and 
Finch  (i®) 

Rasor  and 
McClelland (11) 

Lowenthal 13^ 

Kraftmahker 

2000 

0.1658 

0.1713 

0. 1746 

0.1710 

0.1660 

2200 

0.1683 

0.1754 

0.1807 

0.1764 

0.1692 

2400 

0.1709 

0.1808 

0. 1884 

0.1825 

0. 1723 

2600 

0.1732 

0.1896 

0. 1979 

- 

0.1754 

2800 

0.1756 

0.2026 

0.2111 

- 

0.1786 

3000 

- 

0.2205 

0.2328 

- 

- 

The  dependence  of  the  specific  heat  on  temperature  as  determined  from  enthalpy 
data  is  critically  dependent  on  the  form  of  the  polynomial  chosen  to  represent  the 
enthalpy  data.  Also,  the  results  from  drop  calorimetry  are  subject  to  question 
since  the  reference  state  may  not  be  re- attained  during  the  quenching  operation. 
Pulse  heating  experiments  require  the  accurate  determination  of  the  rate  of  tem¬ 
perature  change  during  transients  but  with  recent  improvements  in  instrumenta¬ 
tion  this  need  not  be  a  major  obstacle.  Thus  pulse  heating  techniques  are  superior 
to  drop  calorimeter  techniques  for  high  temperature  specific  heat  determinations 
for  electrical  .  onductors . 

The  direct  heating  experiment  is  clearly  adaptable  to  specific  heat 
determinations  using  pulse  heating  or  cooling.  The  governing  equation  then  becomes 
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Eq.  (2) 


d 

dZ 


p€„a 

— F—  (T4  -  T  4) 
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L  dT 
^  A  dZ 


6  T 


where  6  is  the  density  and  T  is  the  rate  of  temperature  change.  For  long  rods 
dT/dZ  =  0  and  Eq.  (2)  becomes 


lie 

A2  A 


(T4  -  T  4) 
o 


-  C  5T 
P 


Eq.  (3) 


Previous  users  of  this  general  method  at  high  temperatures  have  used  wire 
samples  thin  tubes  or  an  auxiliary  external  heater  The  use  of 
an  auxiliary  heater  permits  relatively  small  transients  from  an  arbitrary  base 
temperature.  However,  the  use  of  the  heater  greatly  increases  the  apparatus 
complexity  and  the  radiation  loss  problem.  In  this  case  Eq.  (3)  becomes 


I» £ 

A!  A 


(T4  -  T  4)  -  C 


H 


6  T 


Eq.  (4) 


where  a  is  the  total  hemispherical  absorptance  of  the  sample  and  the  term 
aTo4/€H  must  be  determined  accurately.  When  a  heater  is  not  used,  pulse  heating 
to  high  temperature  may  use  large  heating  rates  and  Eq.  (3)  becomes 

^=C  6T  Eq.  (5) 

A2  P 

when  the  radiation  losses  are  made  negligible.  However,  rapid  pulsing  from  room 
temperature  causes  large  temperature  transients  and  possible  changes  in  the 
sample.  Therefore,  for  use  with  the  multiproperty  apparatus  it  was  decided  to 
use  incremental  changes  in  the  power  input  to  achieve  transients.  This  approach 
does  not  require  significant  modifications  to  the  experimental  apparatus  except 
for  the  data  collection  system  and  does  not  put  undue  stress  upon  the  sample  nor 
does  it  involve  quenching. 

One  may,  of  course,  increase  or  decrease  the  rate  of  power  input  and  thus 
observe  either  coding  or  heating  of  the  specimen.  Rates  of  cooling  that  can  be 
achieved  from  a  1/8  inch  diameter  specimen  of  tantalum  by  incrementally  de¬ 
creasing  the  current  are  plotted  in  Figure  4.  From  this  figure  it  can  be  seen  that 
initial  cooling  rates  at  2000  K  are  15  C/sec  for  a  10  percent  change  in  current 
and  increase  to  93  C/sec  for  a  100  percent  decrease  in  power.  These  rates  are 
very  modest  compared  to  the  103  to  105  C/sec  used  by  the  researchers  using 
pulse  heating  at  high  temperatures  without  an  auxiliary  heater  and  should  not 
cause  changes  in  the  sample  structure. 
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III.  EXPERIMENTAL  APPARATUS  AND  PROCEDURES 


The  basic  apparatus  for  steady-state  measurements  is  described  in  the 
literature  5’ 8^.  Briefly  it  consists  of  a  high  vacuum  (10"8  torr)  chamber 
equipped  with  an  optical  window,  regulated  power  supplies,  and  a  guarded  six- 
dial  potentiometer.  The  sample  is  clamped  between  movable  electrodes  inside 
the  vacuum  chamber.  Temperature  measurements  are  made  using  an  automatic 
optical  pyrometer  mounted  on  an  external  platform  which  can  be  raised  or  lowered 
to  view  the  entire  specimen  through  the  optical  window. 

An  overall  view  of  the  apparatus  is  shown  in  Figure  5.  The  automatic 
high  vacuum  system  with  bell  jar  equipped  with  two  rectangular  optical  windows  is 
shown  to  the  left  of  center.  The  automatic  optical  pyrometer  mounted  on  the  plat¬ 
form  is  located  in  front  of  one  of  these  optical  windows.  The  regulated  DC  power 
supplies  are  located  in  panels  near  the  right  side  of  the  picture.  The  guarded  six- 
dial  potentiometer  facility  is  not  visible  in  this  view. 

In  order  to  perform  the  transient  measurements  in  which  the  rate  of  power 
is  incrementally  changed,  it  was  necessary  to  provide  a  means  of  rapidly  changing 
the  rate  of  power  input  and  to  provide  a  means  of  measuring  the  rate  of  tempera¬ 
ture  change.  The  output  of  the  power  supplied  was  changed  by  switching  a  variable 
resistor  in  parallel  with  the  control  resistor.  By  changing  the  value  of  the 
variable  resistor  it  was  possible  to  change  the  current  level  over  a  wide  range  and 
thus  control  the  rate  of  cooling  or  heating.  This  switching  network  is  visible  on 
top  of  one  of  the  regulated  power  supplies  (Figure  5). 

It  was  found  that  the  Kepco  power  supplies  on  hand  could  be  decreased  from 
one  current  setting  to  another  within  250  milliseconds  including  overshoot  when 
the  power  supplies  were  operating  in  the  constant  current  mode.  The  level  would 
change  to  l/e  of  its  final  value  within  100  milliseconds.  The  time  response  was 
about  half  of  this  value  when  the  current  level  was  increased.  The  typical  re¬ 
sponses  are  shown  in  Figure  6.  Thus  by  neglecting  the  data  obtained  during  the 
first  0. 3  seconds,  the  use  of  Kepco  power  supplies  was  feasible. 

Unfortunately,  the  largest  Kepco  supply  available  is  capable  of  delivering 
only  100  amps.  By  paralleling  two  one-hundred  nmp  supplies,  it  was  possible  to 
obtain  200  amps  at  8  volts.  Using  this  arrangement  it  was  possible  to  increment- 
ally  change  the  temperature  of  the  usual  size  metal  samples  (0. 125  in.  dia.  by 
12  in.  long  rods)  up  to  2500  K.  By  paralleling  a  third  Kepco  supply  (on  hard)  and 
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FIGURE  4.  COOLING  RATES  RESULTING  FROM  INCREMENTAL  CHANGE  IN  CURRENT 
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FIGURE  6.  RESPONSE  OF  POWER  SUPPLIES  TO  STEP  CHANGES 


decreasing  the  effective  sample  length  to  10  inches,  it  should  be  possible  to 
attain  3000  K.  Larger  electrical  conducting  cable  is  required  in  order  to  in¬ 
corporate  this  modification. 

The  automatic  optical  pyrometer  used  can  be  operated  in  the  transient 
mode  by  operating  directly  off  the  photomultiplier  tube.  However,  this  signal- 
to-noise  ratio  is  unsatisfactory  and  the  signal  must  be  filtered.  The  response 
time  is  controlled  by  the  filtering  network.  Consequently  a  combination  DC  am¬ 
plifier  and  AC  filter  was  designed  and  constructed  specifically  for  this  transient 
work.  This  electronic  device  is  located  just  below  the  automatic  pyrometer 
(Figure  5).  The  output  of  the  filter  is  fed  into  an  oscilloscope  (shown  in  Figure  5) 
or  into  an  X-Y  recorder  (not  visible  in  Figure  5).  The  sensitivity  of  the  pyro¬ 
meter-amplifier-filter  circuit  is  controlled  by  varying  the  high  voltage  supplied 
to  the  pyrometer  photomultiplier  and  by  varying  the  resistance  in  the  amplifier 
circuit.  Sensitivities  from  5  to  50  mV/C  were  commonly  used.  Several 
neutral-density  filters  were  mounted  on  a  rotating  wheel  on  the  front  of  the  pyrom¬ 
eter  (Figure  5)  and  were  used  while  adjusting  the  sensitivity  of  the  electronic 
output. 


After  attaining  a  satisfactory  sensistivity,  the  electronic  system  was 
calibrated.  This  was  done  either  using  a  recorder  or  an  oscilloscope.  First  the 
transient  output  was  adjusted  to  approximately  the  top  of  the  scale  on  the  record¬ 
ing  device  and  a  trace  of  this  output  was  recorded.  Then  the  current  to  the 
sample,  voltage  drop  across  the  voltage  probes  and  the  pyrometer  output  (ref¬ 
erence  mode)  were  recorded.  The  pyrometer  was  switched  back  to  the  transient 
mode  and  the  sample  current  was  changed  until  the  transient  output  approximated 
one  scale  division  less  on  the  recording  device.  A  trace  was  made  of  this  output 
and  the  current,  voltage  and  pyrometer  output  (reference  mode)  recorded.  This 
procedure  was  repeated  until  the  last  available  scale  mark  on  the  recording  device 
was  used.  A  typical  calibration  is  shown  in  Table  III.  The  sample  current  was 


then  returned  to  its  original  value  ar*d  the  final  current  (Ij)  adjusted  to  a 

desired  value  by  temporarily  switching  in  a  variable  resistor.  After  the  sample 
reached  equilibrium  at  current  was  switched  to  1^  and  the  tempera¬ 

ture  was  recorded  as  a  function  of  time.  Using  the  same  calibration,  other  values 


of  1^  were  chosen  and  additional  traces  were  obtained.  An  oscilloscope  recording 
is  illustrated  in  Figure  7  for  the  calibration  of  Table  III. 
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TABLE  IU 


TYPICAL  CALIBRATION  OF  ELECTRONIC  OUTPUT 
(Tantalum  Sample) 


Calibration 

Mark 

Sample 

Current  (Amps) 

Voltage 
Droptf Volts) 

Pyrometer 
Output  (MV) 

Sample 
Temp.  (K)j! 

1 

117.705 

0.30166 

299.05 

1906 

2 

114.960 

0.29183 

295.45 

1885 

3 

111.947 

0.28121 

291.11 

1860 

4 

108.755 

0.27009 

286.83 

1837 

5 

104.923 

0.25688 

281.61 

1808 

6 

100.320 

0.24127 

275.46 

1774 

^  L  =  2.4623  cm,  A  =  0.072923  cm2,  §  =  16.595  gm/cm3 

>Jc 

From  pyrometer  output  and  emittance  and  window  corrections. 


Using  the  data  of  Table  ID  and  Figure  7  it  is  possible  to  calculate  specific 
heat  values.  The  Table  IB  data  can  be  considered  to  be  a  series  of  equilibrium  values 
for  an  infinitely  long  sample.  Under  these  conditions 

I2  P 

=  P€u  a  (T  4  -  T  4)  Eq.  (6) 

^2  H®  ®  O  ' 

The  data  from  Figure  7  can  be  considered  to  be  for  the  nonequilibrium  infinitely- 
long  sample  (Eq.  3).  At  any  particular  instant  during  the  transient  Eq.  (3)  may  be 
written  as 


6  T 


Substituting  Eq.  (6)  into  Eq.  (7) 


I  2p 

00  ~00 


=  C  6(f) 
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00 


Eq.  (7) 


Eq.  (8) 
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TIME  (SEC.) 

FIGURE  7  TYPICAL  TEMPERATURE  VERSUS  TIME  RECORDING 
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CURRENT  (AMPS) 


Since  at  each  instant  during  the  transient,  0  =  0^  and  T  - 


C 

P 


o  <V  -  O 

A2  6(T)b> 


Eq.  (9) 


From  Table  in  and  the  sample  geometry,  p  may  be  calculated 


_  A  E  A 

L 

00 


Eq.  (10) 


The  calculations  for  the  data  of  Figure  7  and  Table  in  are  given  in  Table  IV. 


All  factors  except  (T)^  entering  into  Eq.  (9)  are  known  to  about  5  significant 
figures  (Table  IV).  Thus  the  accuracy  of  the  specific  heat  determination  is  strictly 
controlled  by  the  accuracy  obtainable  for  (T)^.  Values  for  (T)m  are  obtained  from 
the  slope  of  the  curve  at  various  values  of  each  T^  (Figure  7).  Using  mechanical 
recorders  or  oscilloscopes,  the  usual  accuracy  is  about  5  percent.  This  may  be 
improved  by  about  a  factor  of  two  by  very  careful  and  arduous  efforts.  However,  an 
improvement  by  an  order  of  magnitude  could  be  achieved  using  the  procedures  out¬ 
lined  in  the  following  section. 


TABLE  IV 


SPECIFIC  HEAT  CALCULATION  FOR  DATA  OF  TABLE  HI  AND  FIGURE  7 


Temp.  (K) 

O  (Ohm-cm) 

I J  -  Ij2  (Amps2) 

T  (T/sec) 

Cp*  (J/gmK) 

1860 

74.395  x  10-6 

8051.2 

40.6 

0.167 

1837 

73.  550  x  10~® 

7346.7 

35.5 

0.172 

1808 

72.608  x  10"8 

6527.9 

31.5 

0.170 

p 

A2  6  T 


where  A2  6  =  0.  088238  gm/cm. 
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IV.  EXPERIMENTAL  RESULTS  AND  DISCUSSION 


Experimental  measurements  were  made  on  tantalum  and  L  nt alum-10  tungsten 
alloy.  These  materials  were  chosen  primarily  because  samples  which  were  used 
in  the  apparatus  for  steady  state  measurements  were  on  hand.  Thus,  thermal  con¬ 
ductivity,  electrical  resistivity,  Lorenz  riitio,  total  hemispherical  emittance  at 
0.65  microns  and  Thomson  coefficient  data  were  already  available  Therefore, 
developing  the  capability  of  measuring  specific  heat  on  these  samples  automatically 
demonstrates  the  feasibility  of  calculating  thermal  diffusivity  as  well  as  enthalpy 
from  the  data  and  increases  the  number  of  properties  measured  with  the  multiproperty 
apparatus  to  nine.  In  addition  to  the  convenience  of  having  specimens  available  with 
a  suitable  geometry,  the  specific  heat  of  tantalum  at  moderate  temperatures  (1000  to 
2000  K)  is  relatively  well  established  and  the  specific  heat  of  the  Ta-lOW  alloy  should 
be  very  close  to  that  of  pure  tantalum. 

The  tantalum  samples  are  described  elsewhere  They  consist  of  1/8  inch 
outside  diameter  by  0.040  inch  inside  diameter  thick  wall  tubes  of  metallurgical  grade 
material.  The  nominal  composition  is  99.  9  percent  tantalum  with  small  amounts  of 
silicon,  magnesium,  copper,  and  calcium  and  traces  of  iron,  nickel,  and  zirconium 
present.  The  density  is  16.593  gm/cm3.  The  Ta-.'  V  alloy  is  from  a  batch  used  in 
a  high  temperature  reference  program.  Only  a  few  measurements  were  made  on 
this  material. 

The  original  measurements  were  made  with  the  pyrometer  sighting  on  the 
sample  surface.  It  was  discovered  that  such  a  procedure  yielded  poor  specific 
heat  results  when  data  taken  near  the  start  of  the  transient  were  used.  Similar  re¬ 
sults  were  obtained  with  the  Ta-lOW  alloy  sample.  There  may  be  several  reasons 
for  this  behavior.  The  rate  of  change  of  surface  temperature  just  after  the  initia¬ 
tion  of  the  transient  may  not  be  equal  to  the  rate  of  change  of  the  effective  sample 
temperature  for  a  1/8  inch  or  larger  diameter  rod  of  tantalum.  Also  the  change  in 
power  input  is  not  a  step  function  (Figure  6).  While  a  portion  of  the  data  taken  im¬ 
mediately  after  the  initiation  of  the  transient  was  discarded,  it  is  possible  that  the 
first  part  of  the  data  used  in  the  calculations  were  still  influenced  by  the  departure 
from  a  step  function. 

Some  measurements  were  made  using  the  change  in  surface  temperature 
measured  during  the  middle  and  near  the  end  of  the  transient.  Specific  heat  values 
calculated  from  these  data  were  markedly  improved  and  some  of  these  data  are  in¬ 
cluded  in  Figure  8  which  compares  the  present  results  to  literature  values.  The 
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improvement  in  these  results  is  believed  to  be  caused  primarily  by  rate  of  change 
of  surface  temperature  approaching  the  rate  of  change  of  effective  temperature  near 
the  end  of  the  transient. 

A  0. 060  inch  diameter  hole  was  drilled  perpendicular  to  the  sample  axis  to 
intercept  the  inner  bore.  The  effective  emittance  of  this  hole  was  determined  from 
the  brightness  temperature  and  from  the  electrical  resistivity  versus  temperature 
curve.  Thus  temperature  measurements  of  the  interior  of  the  sample  during  tem¬ 
perature  transients  could  be  directly  measured  using  the  automatic  optical  pyrometer. 
Using  these  measurements,  a  significant  improvement  in  the  specific  heat  results 
was  noted.  Results  calculated  from  measurements  near  the  initiation  and  near  the 
end  of  the  transient  were  in  much  closer  agreement  than  those  obtained  from  surface 
measurements . 

Typical  results  using  temperature  measurements  in  the  hole  are  shown  in 
Figure  8.  From  this  figure  it  is  seen  that  the  present  experimental  results  using 
the  hole  are  in  reasonable  agreement  with  literature  values  ^10,  14, 15\  but 

that  the  scatter  is  considerable  -  up  to  ±  8  percent.  Cooling  rates  from  31  to  90  K/sec 
were  used  in  obtaining  these  results.  Most  of  the  scatter  in  experimental  values  is 
due  to  the  uncertainty  in  determining  T  (Eq.  9).  This  scatter  would  be  greatly 
reduced  with  improved  data  recording  instrumentation. 

Specifically,  if  a  digital  data  acquisition  system  including  a  small  computer 
for  transient  data  storage  (or  an  auxiliary  storage  unit)  were  used,  the  voltage, 
current,  temperature,  and  rate  of  change  of  resistivity  could  be  recorded  during 
the  transient  to  a  much  greater  degree  of  accuracy  than  the  present  system  based 
on  an  oscilloscope.  Then  T  could  be  determined  from  the  rate  of  change  of  resis¬ 
tivity  (p)  and  the  known  temperature  coefficient  of  resistivity,  i.e. ,  T  =  p/(dp/dT). 
This  procedure,  which  was  successfully  used  by  Taylor  and  Finch  has  the 
advantage  of  determining  the  rate  of  change  of  the  effective  temperature,  rather 
than  the  rate  of  change  of  the  surface  temperature.  Also  problems  associated  with 
the  temperature  dependence  of  the  spectral  emittance  and  window  correction  and  the 
non-linear  pyrometer  response  are  avoided.  In  addition,  the  measurement  of  p  can 
be  obtained  quite  accurately  using  a  digital  data  acquisition  system.  For  the  ex¬ 
ample  given  in  Table  III  and  Figure  7,  the  change  in  voltage  at  66.  940  amps  for  a 
temperature  change  from  1906  to  1774  K  is  0.010569  volts.  Using  a  digital  data 
system  and  an  amplifier,  this  voltage  change  could  be  measured  to  better  than 
0.1  percent  and  T  could  be  determined  within  1  percent.  Thus  the  specific  heat 
could  be  determined  within  ±  1  percent  at  high  temperatures. 
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V.  SUMMARY  AND  CONCLUSIONS 


The  multiproperty  apparatus  can  be  used  to  determine  the  specific  heat  at 
high  temperatures  (1000-2800  K)  on  the  same  sample  used  for  thermal  conductivity, 
electrical  resistivity,  total  hemispherical  emittance,  spectral  emittance  at  0.65 
microns,  Thomson  coefficient,  and  Lorenz  ratio  measurements.  Enthalpy  and  ther¬ 
mal  diffusivity  can  thus  be  calculated  from  the  results.  The  specific  heat  deter¬ 
mination  involves  small  transients  from  an  established  steady-state  temperature, 
and  should  not  put  severe  strain  on  the  samples  nor  cause  significant  internal 
changes. 


The  equation  used  to  calculate  specific  heat  from  the  experimental  data  is 


o  <V  -  ■’> 

A*  6  T 


All  quantities  entering  into  this  equation  except  T  could  be  measured  to  about  five 
significant  figures.  Consequently,  the  accuracy  of  the  Cp  determination  was  con¬ 
trolled  by  the  ability  to  determine  T.  Using  an  automatic  pyrometer  and  an  oscil¬ 
loscope,  this  uncertainty  was  about  ±8  percent.  However,  using  a  commerically 
available  digital  data  acquisition  system,  the  inaccuracy  should  be  reduced  by  about 
an  order  of  magnitude. 
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